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Smooth interpolation of lattice gauge fields by signal processing methods 
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We digitally filter the Fourier modes of the link angles of an abelian lattice gauge field which produces the 
Fourier modes of a continuum A l _ l (x) that exactly reproduces the lattice links through their definition as phases of 
finite parallel transport. The constructed interpolation is smooth (C°°), free from transition functions, and gauge 
equivariant. After discussing some properties of this interpolation, we discuss the non-abelian generalization of 
the method, arriving for SU(2), at a Cayley parametrization of the links in terms of the Fourier modes of A^(x). 
We then discuss the use of a maximum entropy type method to address gauge invariance in the non-abelian case. 
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1. Lattices and signal filters 

The explicit relationship between an array of 
links C/ M and the continuum field A^(x) to which 
it corresponds is an important ingredient in com- 
paring lattice results with perturbation theory. 
While the mapping of continuum interpolations 
to lattice configurations is (very) many-to-one, 
there are various simple algorithms which give 
unique continuum fields that serve our purposes 
to lesser or greater degree. Probably the sim- 
plest of these is the piecewise constant interpo- 



lation where links are U u 



a iaA, 



, and A^(x) is 



constant along links. Despite its simplistic prop- 
erties, this interpolation is a mainstay of lattice 
gauge theory. 

The purpose of this contribution is to explore 
better interpolations of a lattice configuration us- 
ing some ideas borrowed from signal processing. 
The starting point is that the interpolating field 
be smooth (C°°), and that it reproduces the lat- 
tice links exactly through the definition of parallel 
transport between x and x + a 



Ufj,(x, x + a) = Ve 



f x+ St A^x 1 ,x 2 ^..J,...^ ) 
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In light of developments in improved actions over 
the past few years, we know that a better rep- 
resentation of the underlying continuum field is 
extremely powerful. Knowing A^(x) allows us 
to compute (classically) "perfect" operators di- 
rectly in the continuum such as F?„, F^F^, 
d^A^, and Dct[p (A^)], providing a greatly im- 



proved action, topological charge density, gauge 
fixing functional, and fermion effective action. It 
also provides a method of implimenting chiral 
fermions del Desperados [HQ- 

The interpolation of an unknown function from 
a discrete set of measurments is the bread and 
butter of digital signal processing j| and I find it 
useful to think in terms of the following analogy 
between the lattice and a digital filter: 



input 



lattice 



filter 



output. 



Thus, a lattice is nothing but a band limited dig- 
ital filter of continuum gauge fields. Ususally a 
signal filter is represented by its transfer function 
Xk in fc-space: 



/in(fc) 



Afc/in(£;) 



/out(fc). 



Similarly the effect of a lattice on A^(x) has a 
simple representation, as follows. 

2. An abelian extrapolation 

An abelian link U^(x : x + a) = e %B is defined 
by an angle — 7r < O^ix + |) < 7r, assigned to 
the point (a; 1 , a; 2 , ...,x^ + ...,x D ). The lattice 
transfer function Xk then follows imediately. We 
first Fourier transform the angle field 
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_ y ' g (fc)e i ' JlxI ~' l"* ! n( :r +f)H hifcoa: 

fc 

= J2 [ X a dt»A^(ky kixl+ - +k ^ + - +ikDxD 

fc 

= ^ \{k^)A l ,(k)e tklxl+ - +k ^ x+ ^ + - +ikDxD 
2 sm(ak fi /2) 



ft 
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This filter is so common in signal analysis that it 
has a special name, sinc(fc); it arises when the 
data are integrals of the underlying signal, rather 
than samples. A familiar example is a CCD cam- 
era where pixels represent the integrated count of 
photons over some length of time. 

We then have our extrapolation: the Fourier 
modes of the continuum field corresponding to 
the lattice with angles 6* M , are 



•MM 



(3) 



2.1. Properties 

• Gauge equivariance 

It is very easy to show that the above con- 
traction is gauge invariant when the continuum 
gauge section G(x) is extracted likewise from its 
lattice counterpart by sampling (ie. \k = 1). 
This gives us a continuum lu(x) whose exponen- 
tial G{xi) — e tu ( Xi ^ at lattice sites Xi is the gauge 
transformation applied to the lattice field. Then 
the following diagram holds: 

Lattice Continuum 
U„(xi + § ) — ► A^x) 



u(x) (4) 



1 



A^{x) - d^u(x) 



G{ Xi ) — > 
1 

+ — A'^x) 

ie. the extrapolation of the gauge transformed 
lattice field V , is identical to the continuum 
gauge transformation A'^ (using the extrapolated 
section oj{x)) of the original extrapolation A^x). 
This is true at arbitrary x. 

• Topological charge 

Since the extrapolated gauge field is C°° and 
periodic, on a manifold without boundaries it 



must have topological charge Q = 0, a sim- 
ple fact maintaining Stokes law. While some 
might be unhappy that an extrapolation has 
Q = when other lattice definitions give nonzero 
charge, I view this as a feature and not a bug. 
We should remember that nontrivial topological 
charge stems from nontrivial transition functions; 
if our lattices are periodic, then our gauge fields 
should have Q = 0. Any non-trivial topologi- 
cal charge arises from a misrepresentation of the 
geometric relationship between A^ and F^, and 
means that we have allowed singular gauge sec- 
tions. 

The topological charge density on the other 
hand, will not be zero on the extrapolation, and 
should be a perfectly good measure of the true 
topological fluctuations and character of the field. 
In particular, FF computed from the extrap- 
olated A^(x) can be integrated over many re- 
gions small compared to the volume of the lattice, 
whose average will then give an accurate measure 
of Q/Vol. 

• Large gauge transformations 

In the diagram (Q), extrapolated gauge trans- 
formations e lu ^ are restricted to those with zero 
winding if we use a strictly periodic Fourier basis. 
However two configurations t/ M and UL(x) which 
are related by a large gauge transformation will 
only differ in their zero momentum modes, and 
the extrapolation to A^ (x) and A'^ (x) will repro- 
duce this difference correctly, despite the fact that 
it cannot produce the continuum gauge transfor- 
mation between them. 

3. Nonabelian extrapolation 

There are two essential difficulties in applying 
the idea of the last section to nonabelian gauge 
fields. 

3.1. Nonabelian parallel transport 

...is very nonlinear. The relationship between 
link and gauge field is given by the path ordered 
exponential 



Ufj,(x, x + a) = Ve 



x+a 
f dt A„ 

•J X 
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rx-\-a 

= 1 + i dtAf,(t) 

J X 

- - dt 2 dtx\®{t 2 - ti)A p (t2)A M (ti) 

+ e(*i-t 2 )A M (ti)A M (ta)} 
+ ••• (5) 

where 0(f) = (a; > 0)? 1 : is a step function, 
and the permutations over ordering are explicitly 
displayed. We must evaluate the Fourier trans- 
form of each term, then do the sum. The main 
part of the computation is evaluating integrals of 
the form 

nx-\-a px-\-a px-\-a 

dt n ... dt 2 dhe(t n - tn-i) . . . 0(ta - h) 

J x J X J X 

(6) 



xe 



ik n t n -\-ik n 



-iH Hkiti 



I have evaluated terms up to third order in A p 
using 



e(t) 



-IT 1 

2ir L ^ in 



e tt +t + 7r 
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on a gauge field with no zero momentum mode. 
There are interesting cancellations which give 
hope that the entire sum can be done analyti- 
cally. For example in the second order term 



h = 



x-\-a nx-\-a 

dta / dhQ(t 2 - t 1 )e ik2t2+ikltl 



A/^+nA/^ — n 



At A 



fci 



/ d d 
\<9fc 2 dk\. 



(8) 



the sum over n ^ (done via the C function con- 
tour) exactly cancels the partial derivative term, 
leaving 



h = ^Afe^+^+t) 



(9) 



which is nicely factorized in k\ and k 2 , linear in 
each Afe, and centered at x+ §. 

The third order term with two nested theta 
functions ©(£3 — £2)0(^2 — ti) is tedious but 



straightforward. While many terms cancel, cer- 
tain derivatives of A& remain in an individual in- 
tegral of type d|) , yielding 



h = 



4Z2 {^fcs^fca^fci ~^fe3Afe 2 A' fel 



7i" 2 Afc 3 Afc 2 A fel je 



i(k 3 +k 2 +k 1 )(x+^ ) 
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However the derivative terms do cancel when 
summed over the 3! permutations of the ordering 
of A ll (k n ) leaving only the last term in (|l0|); again 
factorized in k n , linear in each A&, and centered 
at x ~\~ 2 • 

For SU(2) fields, I have summed the series in 
based on the above results and using explicit 
formulae for the permutations of products of oi\ 
the result is interesting but incomplete. Assum- 
ing that permutations cancel the derivative terms 
at each order, the coefficient of the last term 
Afc„Afc n _ 1 . . . Afcj is easily obtained and produces 
a series which evaluates to the unitary matrix 



1 + ^AfcA^fc) 



1-^AfeA^fc) 

k 
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This is a Cayley representation of U^, similar 
to that recently espoused by PeriwalQ. This 
shows that there must be more terms appear- 
ing at higher orders, since the correct result must 
faithfully represent 



U^a, 6)17^(6, c) = C/ M (a,c), 



(12) 



which (jllj) unfortunately does not do. Nonethe- 
less, it seems an interesting result. 

3.2. Gauge invariance and k c 

In the abelian case we can content ourselves 
with a frequency cutoff in A^(k) and 0{k), and 
due to the linearity of gauge rotations, the entire 
construction is equivariant. In the non- abelian 
case gauge transformations, which involve prod- 
ucts of fields G'{x)Af 1 {x)G{x) — . . ., mix Fourier 
modes of G and A^ beyond the cutoff momentum 
k c . 

Again let's turn to engineering for a solution. 
In situations where one needs to know the power 
spectrum at a frequency other than at discrete 
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values of the measured FFT, the method of Max- 
imum Entropy says that the FFT is actually equal 
to a Laurent series in the complex plane, which 
can be approximated by a Pade polynomial 



/0) = £/(*)** = 



i 



E« 



(13) 



where Zk = e . We can use the Pade approx- 
imation for f(k) at any k : \k\ < k c . 

Another ingredient, known as aliasing, is that 
when doing an FFT, the "measured" values of 
fk include contributions from each multiple of k c . 
Thus, 



(14) 



where F(k) is the measured FFT, and f(k) is the 
actual spectrum. Can we use the aliasing to get 
information about the spectrum at all k? 
Consider the Pade approximation to f(k) 



f(k) 



A(k) 
B(k) 



(15) 



where A and B are polynomials. Including the 
effect of aliasing this becomes (here we use lat- 
tice units a = 1, so that k c = L and the fourier 
frequencies — -| < k < j — 1 are integer). 



F(k) = 



E 



A{k - nL) 
B(k-nL) 



E A{k~zL) 
nCOt{7:z) B(i—^L) 



(16) 
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where the residues are computed at the zeros of 
B(k — zL). If we assume a simple form for B such 
a#| 



B(k) = L 



2L 



■2L 



(17) 



then the poles and residues are easy to find, and 
we have a linear system of equations for the co- 
efficients of A(k). Since A^(x) is real, take A(k) 
real and even for the Re/, and imaginary and odd 



1 Engineers will recognize this as a 2L pole, lowpass But- 
terworth filter. 



for Im/. A(k) is of order k^ 2L 2 ) ensuring finite 
power in f(k). Then, for example 

where 4> n = e l7r< - n+ ^^ L , and A2 P is the coefficient 
of k 2p in A{k), with a similar formula in 2p — 1 
for Imf(k). 

4. Outlook 

I have attempted to outline the issues sur- 
rounding a better interpolation of lattice gauge 
fields. In the nonabelian case there is a large 
technical difficulty in unraveling the nonlinearity 
of parallel transport, which however bears some 
promise of analytic tractability. On the other is- 
sue of gauge equi variance, it would seem we can 
at best make a good approximation to the con- 
tinuous Fourier spectrum of the continuum field 
A fl (x). Notice though that this problem arises 
due to fact that we've used the Fourier basis. 

A major development in signal processing are 
wavelets, which unlike Fourier modes are localized 
both in frequency and space. Perhaps in an ap- 
propriate wavelet basis, the convolution problem 
might not be so bad. In fact the simple piecewise 
linear extrapolation is an expansion in the Haar 
wavelet, one of the simplest examples of a wavelet 
basis. It is gauge invariant precisely because of its 
convolution properties. I hope to explore this ap- 
proach in the future. 

I wish to thank Doug Toussaint, Philippe 
de Forcrand, Yigal Shamir, and Jan Smit for 
thoughtful discussions on the above ideas. This 
work was supported by the US Department of 
Energy, contract DE-FG03-95ER-40906 

REFERENCES 



I. Montvay, |iep-lat/9505015 



W. Bictcnholz and U.-J. Wiese, 



liep-lat/9503022 

S. Stearns, Digital Signal Analysis, Hayden 
Inc. 



4. V. Periwal, tiep-lat/9507007 



